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I.  IIITRODUCTION 


It  is  well  understood  that  tlie  prinary,  if  undeclared,  goal  of  any 
active  radar  systen  is  the  extraction  of  maxinun  target  information. 
Detection  of  the  target  and  estination  of  parameters  such  as  range, 
doppler  and  even  attitude  are  routinely  accomplished.  What  is  not  well 
understood,  however,  is  the  estination  of  a target's  physical  properties 
such  as  size,  shape  and  composition  from  radar  signaling  data.  Tlie 
extraction  of  these  properties  for  naval  vessels  and  their  subsequent 
use  in  identification  or  discrimination  of  the  target  is  the  subject  of 
this  report.  The  identification  of  objects  within  a given  class,  i.e., 
discrimination  and  identification  of  an  object  without  a priori  class 
information,  i.e.,  identification,  are  considered  to  be  separate 
capabilities. 

A basic  understanding  of  the  relationship  between  the  physical 
properties  of  a target  and  its  electromagnetic  scattering  behavior  comes 
from  the  impulse  response  concept  [1].  With  the  normalized  scattered 
field  in  the  far  zone  defined  as 

G(e,((i,p,ja))  = E^(e  ,((>,p,joa)  , (1) 

with  r the  range  and  c the  velocity  of  light,  the  normalized  scattered 
field  and  the  normalized  impulse  response  form  a transform  pair  which 
is 


G(0.'<',P,jw)  = r Fj(e,4i,p,t)e"'^“^dt  . (2) 


Tlie  triplet  e,i}',p  denotes  that  a fixed  target-observer  orientation  and 
fixed  polarization  of  the  incident  field  and  transverse  component  of  the 
scattered  field  have  been  selected.  The  impulse  response  waveform 
"sums  up"  the  scattering  properties  of  an  object  in  one  real  time- 
dependent  waveform.  Thus  specular  contributions,  creeping  waves  and 
other  diffraction  and  reflection  mechanisms  are  clearly  identified. 

Using  convolution,  the  response  to  any  postulated  interrogating  wave- 
form is  obtained.  The  same  approach  yields  the  significant  specular 
range  of  any  particular  scattering  mechanism.  The  role  of  frequency 
becomes  more  evident  from  consideration  of  the  ramp  response  waveform 
defined  as 


G(e.^.P.jw) 

(jl)^ 


j Fp(e,^,p,t)e‘'^“^dt. 


(3) 


1 


r 


Comparing  Equations  (2)  and  (3)  it  is  seen  that  gross  features  of  the 
target  are  defined  at  low  frequencies,  i.e.,  those  frequencies  which 
define  the  ramp  response  waveform,  whereas  successively  more  target 
detail  is  added  by  the  step  and  impulse  response  waveforms.  This  con- 
clusion is  based  on  the  fact  that  for  most  targets 

lim  G(e,(ji,p,ju)  = ACe,4,p),  (4) 

j (*)-»• 

and  for  the  worst  case,  finite  flat  surfaces  viewed  normal  to  the  surface, 
lim  G(J,0,p,jw)  = B{P)ju  . (5) 

Stated  in  the  time  domain,  the  highest  possible  singularity  in  the 
impulse  response  is  an  impulse  (Equation  (4))  for  most  targets  while 
for  a few  targets  (Equation  (5))  a doublet  can  occur.  This  understanding 
of  the  role  played  by  frequency  illustrates  why  modulated  short  pulse- 
type  radar  systems  can  elicite  the  flare  spots  of  a target  but  integrat- 
ing such  information  into  a composite  picture  of  the  target  is,  in  general, 
a very  difficult  task. 

With  any  radar  target  identification  task  there  are  two  basic 
approaches  to  the  problem.  In  the  first  approach  the  starting  point  is 
the  radar  scattering  data  available  from  presently  operational  systems. 
Based  on  such  incomplete  and  nonidentification  dictated  information,  an 
attempt  is  made  to  formulate  methods  and  algorithms  for  identification. 

This  approach  has  not  been  notably  successful.  A second  approach  starts 
with  radar  data  known  to  be  sufficient  for  the  identification  task. 

Note  that  this  does  not  mean  scattering  data  at  all  aspects  and  fre- 
quencies; admitting  however  that  an  exact  inverse  scattering  solution 
probably  exists  with  this  much  information.  In  the  present  context  the 
required  radar  data  are  dictated  by  the  impulse  response  concept.  As 
will  be  seen,  with  this  approach  the  remaining  task  is  to  simplify, 
as  much  as  possible  the  required  complexity  of  an  identification  radar 
system. 

In  this  report  the  required  radar  data  for  conventional  target 
identification  is  reviewed.  Next,  the  special  problems  inherent  to 
naval  vessels  are  enumerated  and  suggested  techniques  for  overcoming 
these  difficulties  are  summarized.  The  suggested  techniques  are  also 
shown  to  be  applicable  to  more  conventional  identification  problems. 

The  next  section  of  the  report  describes  the  reflectivity  ranges  and 
model  targets  used  to  obtain  representative*  radar  scattering  data 

•Representative  in  the  present  context  means  data  from  which  a possibly 
excitation  invariant  set  of  parameters  can  be  obtained.  No  claim  is 
made  (see  section  on  measurement  range)  that  the  data  are  valid  cross 
sections  in  the  usual  sense  of  the  term. 
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for  navel  vessels.  These  data  are  then  presented  together  with  special 
response  waveforms  synthesized  from  the  data.  It  is  also  demonstrated 
that  representative  complex  natural  resonances  can  be  extracted  from  the 
response  waveforms.  A final  section  of  the  report  lists  conclusions 
drawn  from  the  research  to  date  and  recomends  logical  continuations  and 
extensions  of  the  research  effort. 

II.  CONVENTIONAL  TARGET  IDENTIFICATION 

A conventional  method  for  target  identification  [2,3]  is  based  on 
the  fact  that  the  ramp  response  waveform  from  a target  can  be  synthesized 
as 

2N-1 

Fp(6,*,P,t)  “ ^ ^ lG(e,(t.,p,jn(Op)  jcos[nu)^t  + tlnw^)]. 


using  odd  harmonics  or  as 

2N 

’is-  Y. 

° n=l 

using  even  and  odd  harmonics. 
In  Equations  (6)  and  (7) 


|G(e,'J',p,jnu) 


2-  cos[nu)^t  + e(nw^)],  (7) 


5(niu„)  = tan' 


^ Im[G(e,<fr,p,jnWjj)] 
Re[G(e,<j>,p,jnwQ)] 


Thus  if  L is  a maximum  linear  target  dimension  then  discrete  spectral 
scattering  data  for  an  N of  5 and 

..  = 0»47rc 


suffice  to  define  the  ramp  waveform  for  most  targets.  The  result  in 
Equations  (7)  and  (8)  is  coupled  with  the  fact  that  the  physical  optics 
approximation  to  the  ramp  response  waveform  is 

^ ^(2)  , (10) 

TC  c(t-r/c) 
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where  A(z)  is  the  projection  of  the  illuminated  portion  of  the  target 
onto  a plane  perpendicular  to  the  line  of  sight  of  the  radar  (z)  with 
the  illumination  moving  at  1/2  the  free  space  velocity  to  account  for  the 
two  way  path.  It  is  clear  that  the  ramp  response  waveform  contains 
basic  information  on  the  gross  physical  properties  of  the  object.  More- 
over, the  forced  portion  of  the  ramp  response  (Equation  (10))  is  directly 
related  to  the  object's  cross  sectional  area.  This  fact  has  been  used 
to  produce  isometric  three-dimensional  images  of  the  target  [4,5]. 

Conventional  target  discrimination  utilizes  the  fact  that  for  the 
transient  portion  of  the  ramp  response  waveform  a good  approximation  is 

M ^ yAt-t) 

FR(e,.f,p,t)  = ^ Am(e,<(.,p)  e « y(t-t^)  , (11) 

m=l 


where  to  is  the  time  the  interrogating  signal  moves  beyond  the  shadow 
boundary  with  zero  at  the  point  the  response  starts.  Note  in  Equation 
(11)  that  the  excitation  coefficients  contain  the  orientation  and  polari- 
zation dependence  and  the  complex  natural  resonances 


Y 


m 


(12) 


are  only  a function  of  the  size,  shape  and  composition  of  the  target, 
i.e.,  they  are  excitation  invariant.  Methods  for  finding  the  dominant 
complex  natural  resonances  of  a target  from  reaction  integral  equation 
formulations  and  from  transient  signals  or  their  equivalent  spectral 
samples  (Equations  (6)  and  (7))  have  been  discussed  [2, 3, 6, 7].  A 
relatively  new  method  combining  calculated  spectral  data  and  an  integral 
equation  formulation  has  also  been  given  [8].  No  single  method  has 
proven  completely  satisfactory  and  problems  when  noise  and/or  clutter 
contaminate  the  signal  indicate  additional  research  is  needed.  Also, 
while  present  methods  do  not  preclude  use  of  the  complex  natural 
resonances  for  target  discrimination,  they  do  dictate  that  the  discrim- 
ination should  not  be  based  on  extraction  of  complex  natural  resonances 
in  real  time  from  an  unknown  target  signal. 

If  it  is  assumed  that  a library  of  natural  resonances  for  a given 
class  of  objects  is  known,  a priori,  then  discrimination  is  based  on 
prediction-correlation  as 


JT 

n n 


(13) 


L 
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Cr.*l 


where  is  the  measured  ramp  response  and  is  a calculated  ramp 
response.  The  calculated  ramp  response  cones  from  the  general  predictor 
difference  equation  for  transients  [9] 
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N 

f(t)  = ^ tlNAt,  (14) 

n=l 

where  II  is  the  number  of  natural  resonances  and  At  the  sample  interval. 
The  difference  coefficients  in  Equation  04)  are  related  to  the  natural 
resonances  as 


^ ' 


s-e 


y At 
'n 


) = 


n=l 


N 


.N-n 


^N.O  ^ ^ • 


(15) 


thus  given  the  complex  natural  resonances,  yoi  'for  a given  target  the 
difference  coefficients  are  known  for  arbitrary  At.  To  reduce  noise 
problems,  Equation  (14)  is  normalized  by  the  magnitude  of  the  largest 
difference  coefficient.  If  this  is  done,  with  the  N,k  term  say 
having  the  largest  magnitude  then 


F„c(t-k«)  ■ t 


N,n 


n=l 

nA 


iB 


(-1)"^^  FR^(t-nAt), 


N,k 


(16) 


and  each  calculated  point  comes  from,  in  general,  both  prior  and  future 
samples  of  the  measured  waveform.  Numerous  examples  of  successful 
discrimination  using  Equations  (13)  and  06)  have  been  given  [2,3].  In 
the  next  section  certain  general  limitations  and  the  complications 
introduced  by  naval  vessels  are  described.  It  should  be  noted  however 
that  there  is  nothing  magic  about  the  ramp  response  waveform  beyond  the 
discrimination  against  the  influence  of  higher  frequency  scattering 
data.  In  principle,  any  transient  response  waveform  from  the  target 
is  sufficient.  This  assumes  that  a portion  of  the  transient  can  be 
well  approximated  by  a finite  sum  of  exponentials. 

There  have  been  many  applications  of  Prony's  method  and  various 
modifications  of  the  method  reported.  The  references  in  a recent 
dissertation  [10]  are  a good  source  for  Prony's  method.  Two  primary 
and  different  objectives  characterize  the  numerous  applications; 


1 
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1)  good,  finite  exponential  fit  of  time-dependent  signaling 
data,* 

2)  extraction  of  a finite  set  of  complex  natural  resonances. 

In  the  first  application  the  "fit"  to  the  data  is  all  important  - tnere 
is  little  interest  in  the  exponential  arguments.  The  second  application 
stresses  a physical  significance  for  the  exponential  arguments,  i.e., 
only  left  half-plane  arguments  with  negative  definite  real  parts  are 
acceptable.  For  the  identification  problem,  while  the  arguments  or 
complex  natural  resonances  are  a parameter,  it  does  not  necessarily 
follow  that  a physical  significance  must  be  stressed  completely.  Our 
earliest  applications  of  resonances  to  the  identification  problem  used 
ramp  response  estimates  [12].  In  such  data  the  inverse  frequency 
squared  weighting  dictates  only  a few  terms  in  the  exponential  approxi- 
mation for  all  but  very  short  times.  For  these  problems  the  extracted 
resonances  did  have  a physical  meaning.  However  for  complicated 
structures  where  the  responses  are  obtained  from  signal inq  waveforms 
containing  higher  frequencies  (synthetically  or  otherwise)  it  may  not 
be  possible  to  interpret  physically  the  extracted  resonances  nor  may  it 
be  possible  to  insist  on  precise  excitation  invariance.  Thus  it  may  be 
necessary  to  find  that  finite  set  of  complex  numbers  which  best  satisfies 
several  sets  of  data.  This  point  will  be  discussed  later. 

Assume  that  a measured  response  signal  is  obtained  as** 

f^j(t)  = f(t)  + n(t)  , (17) 

where  f(t)  is  the  true  response  and  ji(t)  an  interfering  signal  which 
could  be  noise,  clutter  or  roundoff  error  in  a computer.  The  "measured" 
signal  then  can  be  experimentally  or  computationally  obtained.  Note 
that  this  discussion  is  concerned  with  the  characterization  of  a 
target  and  not  identification,  i.e..  Equation  (17)  are  study  signals 
used  to  obtain  numbers  for  a library.  The  actual  identification 
process  does  not  extract  such  numbers.  The  approximation  is 


fa(t)  = Y.  ^ “f(t)+n(t)  . (18) 

n=l 


*A  frequency  approach  has  also  been  reported,  see  Reference  [11]. 

**Yhe  position,  polarization  and  delay  dependence  is  assumed  to  be 
understood  here  and  hence  is  suppressed. 
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The  difference  equations  inherent  in  the  solution  of  Equation  (18)  are 
M+1  N+1  M+1 


'm  ’ 


1^2  (19) 


m=l  n=l 
where  properly 


m=l 


(20) 


That  is.  Equation  (19)  generates  a system  of  M+1  equations  in  N+1 
unknowns  and  the  error  in  each  equation  is  a function  of  the  sample 
interval,  t2-ti , the  sample  spacing,  a,  the  order  of  the  difference 
equation,  N and  the  number  of  equations  generated,  M+1,  The  Prony  method 
equates  (19)  to  zero  and  sets  C]  to  unity.  In  matrix  form  Equation  (19) 
becomes 


F C = E 


(21) 


where 


i = 1,2, 
j = 1,2, 

•••  M+1 
•••  N+1 

(22) 

^i,j  " ^il 

i = 1,2, 

• • • fl+1 

(23) 

^ij  ~ ^il 

i = 1,2, 

••'M+1  . 

(24) 

and 


In  general  F is  a rectangular  matrix 
M > 2N 


(25) 


with  the  equal  siqn  in  Equation  (25)  yielding  a square  matrix.  The 
inequality  in  Equation  (25)  also  follows  from  Equation  (18)  where  there 
are  2N  unknowns. 

Before  discussing  solutions  of  Equation  (18)  or  equivalently 
Equation  (19)  a difference  in  the  approach  adopted  here  should  be 
noted.  In  Reference  [13]  and  elsewhere  the  exponential  sum  is  based 
on  Singularity  Expansion  Method  (SEM)  ideas.  That  is,  an  impulse 
response  is  written  as 


1 


1 


hj(t)  = 


(26) 


r—'  S t 

L \ “<‘-*1'  ’ 

n=l 


where  it  is  assumed  that  the  singularities  are  simple  poles  (Sp),  the 
transform  of  a possible  entire  function  deleted  and  the  Ap  are  products 
of  the  coupling  coefficients  and  the  natural  modes.  It  is  properly 
recognized  that  the  transient  must  be  delayed  until  the  excitation  is 
removed  (ti>0).  If  this  line  of  reasoning  is  followed  then  for  the 
ramp  response 


h^(t) 


) — e " -e  " ^(1+tS  -t,S  ) 

4 1 j / c \ ^ n I n 


1 (S  ) 
n=l  n' 


>j(t-t^) 


(27) 


Clearly  the  coefficients,  Ap,  can  absorb  the  additional  weighting  in 
Equation  (27)  but  the  additional  terms  in  Equation  (27)  must  also  be 
absorbed  as  exponentials  and  this  in  general  cannot  be  done.  We  have 
for  some  time,  however,  used  the  form  Equation  (26)  to  fit  ramp  response 
waveforms  and  will  continue  the  same  type  of  approximation. 


Returning  to  Equation  (21)  the  solution  of  this  equation  is 
essentially  a departure  point  for  the  two  objectives  mentioned  earlier. 
If  only  a good  fit  to  the  waveform  is  desired  then  Prony's  method  is 
adequate  except  when  signal/noise  (S/N)  power  levels  are  low  [10]*.  If 
Equation  (19)  is  solved  in  the  least  squares  sense  then  the  error  in 
the  fit  to  the  waveform  is  minimized.  This  approach  however  does  not 
optimize  the  exponential  arguments.  Simply  stated,  several  different 
combinations  of  poles  and  residues  can  yield  equally  good  fits  to  the 
response.  The  basic  problem  is  that  the  waveform  fit  is  relatively 
insensitive  to  the  precise  pole  locations. 

Consider  again  Equation  (19),  where  there  are  two  methods  for 
solution  which  have  been  proposed.  The  first  of  these  is  Prony's 
method  where  C]  is  set  to  unity  and  the  remaining  coefficients  are 
chosen  such  that 


M+1 


z 

m=  1 


(28) 


H/e  do  not  Yeel  that  the  S/N  ratio  of  30  dB  (total)  quoted  in  Reference 
[10]  is  realistic  when  the  extraction  process  is  in  non-real  time  and 
not  a one  shot  trial . 
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is  minimized.  This  is  the  standard  least  squares  or  quasi-inversion 
procedure.  A second  method  (eiqenvalue)  is  to  write  Lquation  (19)  in 
the  quadratic  form 


(S)  C . ^ 


where  the  quantity  on  the  right  of  Equation  (29)  is  a minimum. 


/ ■ m ’ 


foi(t+(in+j-l  )A) 


foj(t+(m+i-l  )^), 


are  the  components  of  the  symmetric  matrix  S.  The  coefficients  (C) 
are  then  the  eigenvectors  of  the  minimum  eigenvalue.  It  has  been  shown 
[14]  via  a specific  example  that  the  two  methods  lead  to  different 
squared  errors,  with  the  minimum  obtained  using  the  eigenvalue  approach. 

The  error  parameters  in  Equation  (20)  present  difficulties. 
Unfortunately  these  parameters  are  not  independent,  i.e.,  setting  one 
parameter  imposes  restrictions  on  certain  other  parameters.  For 
example  the  order  of  the  difference  equation,  N,  must  clearly  be  related 
to  the  investigation  interval  t2-ti . If  the  interval  t2-t]  is  pro- 
gressively expanded  to  include  earlier  times  it  follows  that  additional 
exponential  terms  may  be  needed.  Conversely  if  the  difference  equation 
order  is  fixed  and  not  over  specified  then  the  interval  t2-ti  should  be 
bounded.  It  has  been  suggested  [15]  that  the  eigenvalue  approach  can 
be  used  to  establish  N.  This  however  is  based  on  t2-t-|  being  fixed 
and  the  use  of  noise-free  data.  With  noise,  clutter  and  roundoff 
error  present,  as  is  always  the  case,  this  approach  rapidly  deteriorates. 
There  is  also  the  problem  that  if  N becomes  large  discrimination  pro- 
cedures become  unwieldy  [16].  Thus  it  is  advantageous  to  keep  N as 
small  as  possible.  We  postulate  that  N is  best  determined  from  trials 
on  various  data  sets  corresponding  to  the  same  target  which  precludes 
real  time  resonance  extraction  but  is  compatible  with  a predictor- 
correlator  approach. 


The  sample  interval  ^ must  be  carefully  considered, 
sufficient  to  use  the  Nyquist  rate 


It  is  not 


» 


where  Fn,ax  is  the  highest  frequency  in  fo(t).  This  will  yield  an  upper 
bound  but  not  a lower  bound.  Numerous  examples  with  prediction- 
correlation,  where  the  same  difference  equation  as  in  Equation  (19)  is 
used,  have  shown  that  ^ also  has  a lower  bound 
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mm 
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max 


(33) 


If  Fourier  synthesis  is  used  to  obtain  the  transient  then  prior  results 
[16]  indicate 


0.02  1 0.035 

0 


(34) 


but  this  is  based  on  a particular  interval 

t2  - t,  = Tq  - 0.33T^  . (35) 

However,  even  within  the  acceptable  limits.  Equation  (34),  the  results 
are  not  constant,  i.e.,  the  coefficients  in  Equation  (19)  are  a 
function  of  a.  Simply  stated,  the  rectangular  matrix  generated  by 
Equation  (19)  is  poorly  conditioned  numerically.  One  attempt  at  a 
partitioning  scheme  to  improve  the  conditioning  [11]  (a  frequency 
domain  approach)  was  unsuccessful. 

Our  best  results  to  date  have  been  obtained  using  programs  which 
systematically  vary  each  of  the  error  parameters  [16].  When  applied 
to  numerous  data  sets  for  a given  target,  the  results  are  acceptable. 
The  formulation  in  [14]  uses  a Z-domain  approach  and  eigenvalue 
solution,  but  the  Z-domain  at  least  appears  to  be  incidental.  On  the 
present  program  research  has  not  yet  reached  the  stage  where  the 
nuances  of  Prony's  method  are  of  vital  importance.  It  is  sufficient 
that  characteristic  complex  numbers  can  be  extracted  from  the  data  and 
that  these  numbers  are,  via  prediction-correlation,  an  adequate  repre- 
sentation of  the  data.  It  should  be  noted  that  for  the  matched  filter 
responses  of  naval  vessels  the  extracted  (Prony)  resonances  have 
oscillatory  parts  in  excess  of  the  highest  input  frequency.  This  is 
not  surprising  because  the  purpose  of  the  matched  filter  response  is 
to  obtain  a rapidly  damped  response. 

III.  SPECIAL  PROBLEMS  OF  NAVAL  VESSELS 

It  is  evident  from  Equation  (9)  and  the  discussion  on  conventional 
target  discrimination  that  relatively  low  frequencies  are  required 
for  the  ramp  excitation.  The  reason  for  this  is  that  dominant  complex 
natural  resonances  have  been  emphasized.  For  example,  aircraft  target 
natural  resonances  are  used  where  the  wing  and  fuselage  lengths  are 
approximately  one-half  wavelength  at  the  oscillatory  part  of  the 
natural  complex  frequency.  Therefore  the  required  excitation 
frequencies  for  the  transient  response  are  quite  low.  Also,  using 
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Equation  (6)  or  Equation  (7),  complex  radar  scattering  data,  i.e., 
amplitude  and  phase,  are  required.  For  naval  vessels,  if  the  hull 
length  is  considered  the  dominant  linear  dimension,  the  required  fre- 
quencies are  very  low.  For  example  destroyers,  with  a nominal  hull 
length  of  130  meters  would  require  a fundamental  excitation  frequency 
of  0.462  MHz.  Thus  one  special  problem  of  naval  vessels  is  very 
evident. 

Assuming  that  singularities  of  the  target  scattering  function  can 
continue  to  be  modeled  as  simple  poles  at  higher  frequencies,  there  is 
no  basic  reason  why  higher  frequencies  cannot  be  used.  Clearly  as 
frequency  increases  smaller  target  dimensions  approach  an  oscillatory 
resonant  length.  Thus  target  identification  based  on  what  can  be 
termed  "substructure  resonances"  is  a possibility.  A number  of  compli- 
cations can  be  anticipated.  As  already  noted,  the  substructure  reson- 
ances must  be  well  modeled  by  simple  poles.  It  must  also  follow  that 
in  a particular  spectral  range  such  substructure  resonances  must  be 
dominant.  That  is  poles  associated,  for  example,  with  certain  super- 
structure features  must  have  real  parts  at  least  comparable  to  those 
of  larger  structural  features.  Using  an  aircraft  example,  vertical  and 
horizontal  stabilizers  must  have  resonances  similar  to  those  of  the 
wings  and  fuselage  but  at  higher  oscillatory  frequencies.  To  be  useful 
the  damping  part  of  such  resonances  should  not  be  much  greater  than 
those  dominant  for  the  wings  and  fuselage.  It  also  follows  that  for  a 
complex  structure  the  density  of  complex  natural  resonance  locations 
must  increase  as  the  oscillatory  frequency  is  increased.  That  is, 
multiple  resonances  of  longer  structures  may  be  in  close  proximity  to 
the  first  or  dominant  resonance  of  a smaller  structure.  For  this  reason 
it  may  be  necessary  to  use  a larger  number  of  poles  to  describe  the 
transient  response  than  that  necessary  for  the  ramp  response. 
Alternatively,  special  interrogating  signals  based  on  knowledge  of  the 
higher  order  resonances  may  be  needed. 

It  is  well  known  that  radar  signals  from  naval  vessels  exist  in  a 
clutter  background  produced  by  the  sea.  A rare  exception  can  occur 
when  radar  beamwidths  are  such  that  only  portions  of  the  vessel  are 
effectively  illuminated.  An  explanation  of  the  possible  complicated 
electromagnetic  interactions  between  the  vessel  and  the  sea  is  beyond 
the  scope  of  this  study.  In  fact,  a valid  mathematical  model  for  the 
complex  radar  return  from  the  sea  alone  at  arbitrary  frequencies  is  not 
yet  available.  Admittedly  some  success  in  predicting  sea  return  at 
particular  frequencies  has  been  reported,  but  a general  model  has  not 
emerged.  The  sea  then  is  an  added  complication  for  naval  vessels. 

Note  also  that  the  sea  presence  essentially  precludes  obtaining  good 
radar  return  data  from  models  of  naval  vessels  for  higher  full  scale 
frequencies  because  a physical  model  of  the  sea  surface  becomes 
prohibitively  difficult.  For  our  purpose  we  assume  that  if  the  full 
scale  frequencies  are  such  that* 


*The  highest  permissible  frequency  is  somewhat  arbitrary,  but  it  is  felt 
that  Equation  (36)  is  reasonable  if  possibly  optimistic. 
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Fp3  1 25  MHz, 


(36) 


then  the  sea  can  be  modeled  electromagnetically  as  a perfectly  con- 
ducting ground  plane.  With  a highest  model  frequency  of  (see  next 
section) 

= 10.8  GHz  , (37) 

a scale  factor  for  the  model  of  432  or  greater  is  required.  Actually, 
the  choice  of  scale  factors  was  dictated  by  that  available  from 
suppliers  of  plastic  models.*  It  was  decided  that  silver-painted 
plastic  models,  despite  the  size  restrictions,  were  more  realistic  than 
very  crude  homemade  versions  of  the  vessels.  It  is  apparent  however 
that  in  order  to  investigate  full  scale  frequencies  higher  than  25  MHz, 
modeling  would  be  a most  difficult  problem.  The  dilemma  is  serious 
because  both  mathematical  and  experimental  models  of  the  clutter  at  full 
scale  frequencies  greater  than  25  MHz  can  be  seriously  questioned.  We 
have  circumvented  the  problem  for  the  moment  by  restricting  initial 
attention  to  the  lower  frequencies.  The  rationale  being  that  vessel 
identification  should  first  be  established  in  spectral  ranges  where 
reasonably  good  scattering  data  can  be  obtained  on  a model  basis. 

IV.  MULTIPLE  FREQUENCY  REFLECTIVITY  FACILITY 

AND  MODEL  TARGETS 

The  nine  model  targets  used  to  obtain  the  reflectivity  measure- 
ments reported  in  this  study  are  tabulated  in  Table  I.  Only  the  nominal 
lengths  of  the  models  are  given,  other  pertinent  dimensions  can  be 
found  in  Jane's  Ships  of  the  World,  1973.  Admittedly,  minor  differences 
can  and  do  exist  between  the  vessels  reported  in  the  reference  and  the 
actual  models  used  in  this  study.  However,  over  the  span  of  frequencies 
covered  minor  structural  differences  should  not  produce  appreciable 
changes  in  the  scattered  field.  A possibly  more  serious  problem  is 
the  maze  of  antennas  known  to  exist  on  modern  naval  vessels.  Long, 
wire-like  communications  antennas  could  seriously  alter  the  scattered 
field.  No  attempt  has  been  made  to  include  antenna  structures  on  the 
model  targets.  In  this  initial  phase  of  the  study  it  was  felt  that 
principles  could  be  explored  without  the  necessity  of  detailed  antenna 
modeling.  This  amounts  to  assuming  that  antenna  structures  do  not 
dominate  the  scattered  field  over  the  span  of  frequencies  covered. 

There  are  obviously  complex  natural  resonances  associated  with  the 
antennas  but  these  resonances  in  general  would  be  a poor  identification 
tool  because  the  antenna  structures  would  be  approximately  the  same 
size  on  most  vessels.  Associated  vessel  structures  would  alter  the 


*We  are  indebted  to  br.  Arthur  Jordon  of  the  Naval  Research  Laboratory 
for  supplying,  on  a loan  basis,  three  good  500:1  scale  factor  metal 
and  wood  models  of  naval  vessels. 
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resonance  locations  slightly  between  vessel  classes  but  probably  not 
sufficiently  to  affect  identification.  Photographs  of  the  node!  vessels 
on  the  ground  plane  are  shown  in  Figures  1 and  2.  In  every  case,  non- 
metallic  portions  of  the  models  were  silver  painted  to  ensure  realistic 
scattering  data.  As  noted,  the  targets  used  in  this  study  were  for 
the  most  part  commercially  available  scale-model  plastic  ships.  The 
models  were  waterline,  i.e.,  simulating  a floating  structure.  The 
models  were  spray  painted  with  silver-loaded  paint  to  simulate  con- 
ducting metal  targets.  A list  of  the  models  used  and  their  major 
characteristics  is  given  in  Table  I.  While  some  of  these  targets  are 
not  United  States  Na'y  ships  or  even  fairly  recent  foreign  ships,  their 
major  features  are  similar  to  those  of  ships  in  use  today  and  it  was 
felt  that  they  would  be  adequate  for  our  purpose. 

The  measurement  facility  used  to  obtain  data  on  the  ship  models  is 
a unique  bistatic  radar  system  using  the  first  ten  harmonics  of  a 
harmonically  rich  crystal  controlled  1.085  GHz  source,  fo(fo)10fo* 
data  collection  is  performed  under  computer  control  using  an  IBM 
Minimal  Informer.  A simplified  block  diagram  of  the  system  is  shown  in 
Figure  3.  Once  the  data  taking  sequence  is  initiated,  requests  are  made 
to  the  operator  for  the  parameters  of  the  range  geometry.  The  computer 
uses  these  data  to  control  a stepping  motor  which  moves  the  target  one- 
half  wavelength  along  the  range  vector  from  the  target  to  either  antenna. 
Thus  the  two  way  path  differential  is  one  wavelength  at  the  fundamental. 
At  the  fundamental,  this  results  in  one  cycle  of  a slow-doppler  inter- 
ference pattern,  two  cycles  at  the  second  harmonic,  etc.  Since  several 
frequencies  are  used  on  any  one  run  the  return  contains  a superposition 
of  the  interference  patterns.  The  individual  components  are  separated 
by  multiplying  the  sample  values  by  sin(nx)  and  summing  for  n=l,2,---10. 
This  is  also  done  for  cos(nx).  In  essence  this  yields  the  Fourier 
Coefficients  of  a ten  frequency  periodic  pattern.  These  complex  coef- 
ficients therefore  contain  the  magnitude  and  phase  information  on  the 
target. 

For  the  scattering  coefficients  of  an  unknown  target,  three  sets 
of  data  are  taken.  These  are: 

(1)  No  target  run  CNTu^) 

(2)  Unknown  target  run  (U^^) 

(3)  Reference  (conducting  hemisphere)  target  run  (R^). 

where  the  M subscript  denotes  measured  values.  Also,  a calculated  set 
of  scattering  coefficients  for  the  reference  target  at  the  ten  fre- 
quencies is  required  (designated  Rc).  By  the  linearity  property  of  the 
Fourier  Series  the  NTm  term  can  be  subtracted  from  both  Ufi  and  Rfi 
since  it  represents  the  stationary  clutter  return.  By  forming  the  ratio 
Rq  over  Rm  a set  of  complex  calibration  coefficients  is  formed.  The 
calibrated  unknown  target  scattering  coefficients  (designated  Uc)  are 
then  found  by  multiplying  by  the  calibration  coefficients.  This  is 
simply  shown  by: 
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Figure  2.  -1/700  scale  models. 
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Figure  3— Multi  frequency  system  configuration. 
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Using  this  method,  ten  frequency  scattering  data  (both  magnitude  and 
phase)  for  a given  target  and  aspect  can  be  quickly  obtained. 

System  accuracy  is  determined  by  running  another  target  (designated 
the  check  target)  for  which  scattering  data  are  known.  Acceptable  data 
are  usually  based  on  allowable  errors  of  ^5%  in  amplitude  and  ^10%  in 
phase  for  the  check  target.  For  matched  filter  responses  which 
alleviate  the  foldover  errors,  the  10%  phase  error  limit  was  relaxed 
since  phase  data  are  not  used.  For  the  ground  plane  range,  hemispheres 
were  used  for  both  the  reference  and  check  targets.  These  were  modeled 
as  perfectly  conducting  hemispheres  on  an  infinite  perfectly  conducting 
ground  plane.  A far  field  source  and  a far  field  receiver  both  having 
theta  (vertical)  polarization  were  assumed.  The  reference  coordinate 
origin  is  at  the  hemisphere  center  (on  tlie  ground  plane)  with  the  z-axis 
perpendicular  to  the  ground  plane  and  the  x-axis  in  the  plane  bisecting 
the  radar  bistatic  angle.  Tlie  hemisphere  was  then  modeled,  using 
image  theory,  as  a sphere  with  an  image  transmitting  antenna,  see 
Figure  4.  A computer  program  was  written  to  generate  the  required 
magnitude  and  phase  data.  The  magnitude  data  were  output  as  the 
square-root  of  the  radar  cross  section  which  is  proportional  to  scat- 
tered field  strength.  The  phase  data  are  referenced  to  the  center  of 
the  hemisphere  at  its  base.  It  is  for  this  reason  that  all  included 
ramp  response  v^aveforms  for  the  hemispheres  appear  to  start  at  a 
maximum  which  corresponds  to  the  incident  plane  wave  passing  through 
the  phase  center  and  also  intersecting  the  greatest  cross-sectional 
area. 

One  major  problem  with  the  present  system  is  violation  of  the 
far-field  approximation  at  all  but  the  lowest  harmonics.  Usino  a crude 
2 d2/a  approximation  for  the  n-th  harmonic,  the  far-field  range  is 
approximately 


The  antennas  used  are  3 foot  parabolas  with  horn  feeds  (see  Figures  5, 
6 and  7).  With  the  fundamental  of  1.085  GHz  (Xq  = 10.88  inches)  this 
yields : 

r^^  ="  n(19.6  ft)  . (40) 


With  the  present  range  of  about  10  feet  it  is  clear  that  the  targets 
are  not  in  the  far-field.  Wote  that  the  tenth  harmonic  would  require 
a range  greater  than  196  feet  - a very  difficult  requirement.  Our 
purpose,  however,  is  not  to  produce  cross  section  data.  The  ship 
parameters  we  seek  are  independent  of  the  excitation  and  therefore  any 
excitation  which  produces  an  extractable  transient  is  sufficient.  To 
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instrument  a model  range  satisfying  the  above  requirements  would  require 
major  changes  including  TWT  or  equivalent  amplifiers.  This  type  of 
expenditure  was  not  planned  and  is  not  considered  necessary  for  the 
program. 

V.  MHASURED  SCATTERING  DATA  OF  NAVAL  VESSELS 

The  data  obtained  using  tlie  Multi-frequency  System  are  given  in 
magnitude-phase  form.  The  unit  of  magnitude  is  the  square-root  of  cm^ 
yielding  cm.  Thus  the  magnitude  is  proportional  to  the  received 
electric  field  strength.  The  unit  of  phase  is  the  degree.  Phase  is 
measured  from  the  center  of  the  reference  hemisphere  at  its  base.  As 
stated  elsewhere,  the  phase  data  obtained  have  not  been  used  extensively 
due  to  the  use  of  the  matched  filter  step  weighted  responses.  These 
response  waveforms  are  constructed  using  a finite  cosine  series  of  the 
magnitudes  squared  with  a l/n^  weighting  factor  on  the  n-th  term  and  with 
the  phases  set  equal  to  zero.  For  this  reason  the  phase  data  have  not 
been  included  in  this  report,  although  they  are  implicit  in  the  synthetic 
ramp  responses.  The  amplitude  data  are  presented  in  graphical  form  in 
Figures  8 through  37.  For  comparison  an  interpolated  spectrum  is 
given  with  each  set  of  data.  It  has  been  assumed  that  the  spectra 
have  zero  value  at  dc  and  the  eleventh  harmonic.  This  grossly  simulates 
the  lowpass  nature  of  our  finite  cosine  series  and  Rayleigh  scattering 
as  the  frequency  approaches  zero.  The  interpolation  process  is  given 
by 


N = 11  for  ten  harmonics  (41) 


where  G(jnwo)  is  the  complex  scattering  coefficient  at  the  n-th 
harmonic  and  the  Cp  are  modified  impulse  coefficients 
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It  is  not  intended  to  imply  that  the  interpolated  curves  can  necessarily 
be  taken  as  valid  estimates  of  the  scattered  field  amplitude  between 
test  frequencies.  The  data  curves  shown  are  simply  intended  as  a more 
meaningful  display  than  tables  of  amplitude  data  versus  frequency.  The 
curves  obtained  via  Equation  (41)  were  also  found  to  yield  smoother 
results  than  rational  function  fits  of  the  same  amplitude  data  (crosses 
in  the  figures).  A more  complicated  version  of  Equation  (41)  can  also 
be  used  to  estimate  continuous  amplitude  and  phase  spectra  with  the 
added  capability  of  simulating  "chirp-type"  responses.  This  type  of 
processing  is  currently  under  study. 

VI.  PROCESSED  RESPONSE  WAVEFORMS  OF  NAVAL  VESSELS 

The  tabulated  electrical  sizes  of  the  model  ships  given  in  Table  I 
make  it  clear  that  true  ramp  response  waveforms  cannot  be  synthesized 
from  the  measured  scattering  data.  Simply  stated,  the  fundamental 
period  is  such  that  the  transient  signal  does  not  decay  to  negligible 
values  within  this  period  (Fourier  synthesis  automatically  involves 
periodic  waveforms).  Therefore  foldover  errors  occur  in  the  synthesized 
waveform.  The  profile  or  area  function  interpretation  (Equation  (10)) 
cannot  be  made  and  it  does  not  necessarily  follow  that  complex  natural 
resonances  extracted  from  the  waveform  (Equation  (11))  can  be  related 
to  the  physical  target  structure  in  a simple  manner.  It  is  noted 
however,  as  will  be  shown,  that  the  ramp  response  waveforms  obtained 
from  Equation  (6)  or  Equation  (7)  are  characteristic  of  the  target. 

It  is  instructive  to  examine  the  near  foldover  error  for  a calibration 
target;  in  this  case  a conducting  hemisphere  on  a conducting  ground 
plane.  For  convenience,  all  of  the  processed  response  waveforms  given 
in  this  report  are  normalized  such  that  the  peak  magnitude  of  the 
waveform  is  unity.  A vertical  scale  factor  given  on  each  figure  is 
used  to  determine  the  actual  peak  value  of  the  waveform.  The  abscissa 
scale  can  be  given  equivalently  in  either  length  or  time.  For  our 
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purpose  the  scale  is  in  samples  of  the  fundamental  period  (- 


1 


r) 


'1.085x10^ 

and  400  aamples  span  the  fundamental  period.  For  ramp  waveform 
synthesis,  the  ordinate  scale  is  in  units  of  centimeters  normalized  by  n. 
TIte  synthesized  ramp  response  waveform  of  a 1.5  inch  diameter  conducting 
hemisphere  on  a ground  plane  is  shown  in  Figure  38.  The  waveform  was 
synthesized  from  calculated  10  frequency  harmonic  scattering  data 
corresponding  to  the  reflectivity  facility  (fundamental  frequency  of 
1.085  GHz).  Note  that  the  waveform  assumes  a constant  value  over  a 
portion  of  the  fundamental  period  showing  that  foldover  problems  are 
not  present.  Because  the  phase  reference  is  at  the  center  of  the 
sphere,  the  peak  response  happens  to  occur  at  the  start  and  end  of  the 
period.  The  same  waveform,  shifted  such  that  the  ramp  response 
actually  starts  at  the  beginning  of  the  period  and  with  the  constant 
term  removed  is  shown  in  Figure  39.  The  negative  sign  in  Equation  (7) 
has  also  been  included.  It  would  be  a simple  matter  to  calibrate  the 
response  in  Figure  39  in  terms  of  the  profile  function  (Equation  (10)). 
However,  if  the  diameter  of  the  hemisphere  exceeds  the  imposed  limit 
(Equation  (9))  this  interpretation  is  no  longer  possible.  The  near 
maximum  hemisphere  diameter  case  is  illustrated  in  Figure  40,  where 
the  synthesized  ramp  response  waveform  of  a 2.5  inch  diameter  hemisphere 
is  shown.  Figure  40  should  be  compared  witli  Figure  38  and  Figure  41 
with  Figure  39.  The  waveform  in  Figure  40  is  still  a true  ramp 
response  for  a hemispherical  target  on  a ground  plane  but  note  that 
without  the  aid  of  Fiqure  39  the  start  of  the  waveform  is  questionable 
because  a constant  return  over  some  portion  of  the  period  is  not  readily 
obtained.  Increasing  the  diameter  of  the  hemisphere  beyond  2.5  inches 
will  yield  other  characteristic  waveforms,  all  different  in  general 
from  those  in  Figures  39  and  40.  Decreasinn  the  electrical  size  of  the 
hemisphere  below  that  shown  in  Finure  40  will  yield  waveforms  of  the 
same  shape  as  shown  in  Figure  40.  An  exception  to  this  occurs  when  all 
frequency  samples  are  in  the  Rayleigh  region.  In  this  case  the 
synthesized  waveform  will  depart  from  the  anticipated  shape.  In  fact 
a sum  of  equally  weighted  cosine  harmonics  is  obtained.  Finally  it  is 
clear  by  comparing  Figures  38  and  40  that  the  maximum  hemisphere 
diameter  which  will  yield  a true  ramp  response  waveform  is  roughly 


d„,^  < 2.5  inches, 
max  — 


(43) 


The  above  remarks  are  perhaps  obvious  to  one  well-versed  in  Fourier 
synthesis  but  may  need  clarification  for  other  engineers. 

For  further  discussion  of  these  concepts  the  reader  is  referred  to 
any  general  text  on  Fourier  synthesis  (Reference  [17]).  Selected  ramp 
response  waveforms  of  the  model  naval  vessels  tested  are  given  in  this 
section  to  illustrate  that  characteristic  waveform  signatures  are 
obtained.  The  phase  reference  for  the  scattering  data  used  to  synthe- 
size these  waveforms  is  normally  at  the  center  and  base  of  each  hull 
length.  With  this  information  one  can  estimate  the  point  at  which  each 
response  waveform  should  start  using  the  hemisphere  results.  For 
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example,  the  response  of  the  1.5  inch  hemisphere  starts  at  roughly  340 
samples.  The  response  of  the  2.5  inch  diameter  hemisphere  starts  at 
310  samples.  Thus  a length  increase  ratio  of  2. 5/1. 5 yields  a reduction 
of  30  samples  in  the  starting  point.  Estimated  starting  points  are 
indicated  on  the  various  ship  responses,  based  on  the  appropriate  hull 
length  or  beam  width  of  the  models.  If  the  ship  responses  were  true 
ramp  responses  then  a constant  level  could  be  removed  from  the  waveforms. 

In  the  present  case  this  cannot  be  done  and  each  waveform  must  be 
examined  with  an  understanding  that  an  unknown  vertical  shift  is  needed. 

The  starting  points  indicated  for  the  ship  waveforms  are  obviously  gross 
estimates  because  of  the  foldover  errors. 

In  Figures  42  through  49  synthesized  ramp  response  waveforms  of  the 
vessels  given  in  Table  I are  presented  in  the  order  listed  in  Table  I.* 

Each  figure  has  three  waveforms  (a,  b and  c)  corresponding  respectively 
to  bow,  stern  and  starboard  beam  incidence.  These  waveforms  are  included 
primarily  to  illustrate  that  in  most  cases  characteristic  waveforms 
are  obtained.  There  is  no  way  that  a detailed  analysis  of  the  waveforms 
can  be  made.  Ilote,  however,  tliat  in  some  of  the  bow  and  stern  responses 
one  can  discern  an  approach  to  a reasonably  constant  level  from  the 
origin  and  then  larger  magnitude  returns.  This  would  be  anticipated 
as  the  incident  signal  washes  over  the  relatively  clean  hull  and  then 
encounters  the  superstructure.  The  beam  waveforms  differ  in  that  the 
superstructure  is  encountered  immediately. 

From  a complex  natural  resonance  viewpoint  the  ramp  response  wave- 
forms are  not  convenient.  Admittedly  for  most  waveforms  there  is  a 
portion  of  the  return  which  could  be  approximated  by  a finite,  decaying 
exponential  sum.  But  there  is  little  or  no  consistency  among  the 
various  waveforms  and  consequently  less  likelihood  of  obtaining 
excitation  invariant  parameters.  For  this  reason  a different  type  of 
response  waveform  was  desirable.  It  was  also  known  that  the  accuracy  of 
the  recorded  phase  data  was  not  good.  A synthetic  matched  filter- type 
response  can  be  written  as 

2N 

/ lG(e,(ii,p,jnui  ) I |T(e,(j.,p,jn(;L,  )lcos[nw  t+5(nu  )+6(n(jj  )] 

L— t 0 0 0 0 0 

n=l 

(44) 

where  G/C  and  T/6  are  respectively  the  scattered  field  and  filter  in 
phasor  notation.  Discarding  tiie  phase  data  and  using  a step  weighting 


fn^(9.<»'.P,t)  = 


♦The  Iowa  (1:700)  is  not  included  since  it  is  identical  in  construction 
to  the  Missouri. 
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Figure  42(c).  Midway  (1/500)  - Starboard;  Scale:  2.265. 
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Figure  44(b),  Sverdlov  (1/500)  - Stern;  Scale:  0.506. 
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Figure  44(c).  Sverdlov  (1/500)  - Starboard;  Scale; 


1 .665. 
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Figure  45(a).  Missouri  (1/ 
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Figure  45(b).  Missouri  (1/7' 
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Figure  47(c).  Mogami  (1/700)  - Starboard;  Scale: 


900 


0.948. 
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Figure  49(a).  Shimokaze  (1/700)  - Bow  on;  Scale:  0.066. 
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kaze  (1/700)  - Stern;  Scale:  0.079 
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(45) 


H |G(e,(t.,p,jnu)  ) , 

""  2^  2 [nw^t]  . 

n=l  " 

Note  that  the  response  in  Equation  (45)  utilizes  only  the  amplitudes 
squared  of  the  measured  harmonic  scattering  data.  No  claim  is  made 
that  the  matched  filter-type  response  is  an  ideal  identification  tool. 

It  is,  however,  one  method  for  utilizing  widely  separated,  in  frequency, 
measured  scattering  data  whose  phases  are  questionable.  Also,  as  will 
be  seen,  the  approach  yields  response  waveforms  more  nearly  compatible 
with  a complex  natural  resonance  approach  in  the  sense  of  effecting 
decays.  Matched  filter-type  responses  for  the  vessels  in  Table  I are 
shown  in  order  in  Figures  50  through  57  where,  as  before,  three  wave- 
forms are  shown  in  each  figure.  The  filter  responses  are  shown  as 
solid  curves;  the  dashed  curves  are  the  responses  obtained  from  finite 
exponential  approximations  of  the  indicated  portions  of  the  solid 
curves.  We  will  discuss  the  exponential  fits  briefly  at  the  end  of  this 
section,  indicating  that  the  exponentials  obtained  must  be  interpreted 
carefully.  It  can  be  seen  from  the  matched  filter-type  responses 
however  that  good  exponential  fits  can  be  obtained. 

A summary  of  the  complex  exponentials  obtained  from  approximate 
fits  of  the  matched  filter  responses  is  given  in  Table  II,  The 
response  in  Equation  (45)  can  be  viewed  as  an  approximation  to  a 
matched  step  response.  It  must  be  remembered  however  that  calculations 
as  in  Equation  (45)  are  using  a very  few  samples.  It  can  be  anticipated 
that  an  exponential  approximation  of  Equation  (45)  will  require  oscil- 
latory frequencies  larger  than  in  the  original  samples.  It  is  maintained 
that  the  complex  numbers  given  in  Table  II  are  characteristic  of  the 
vessels  but  no  attempt  to  associate  the  numbers  with  particular 
component  structures  is  made.  In  a following  section  it  is  shown  that 
discrimination  of  the  vessels  can  be  achieved  using  the  numbers  in 
Table  II  and  a predictor-correlator  processing.  However,  Table  II 
shows  that  the  numbers  are  not  yet  excitation  invariant.  Swept  frequency 
measurements  over  the  same  bandwidth  covered  by  the  discrete  samples 
would  be  desirable  and  are  planned  for  the  future.  An  example  of  one 
such  measurement  has  been  made.  Two  1/700  scale  models  of  the  Missouri 
were  connected  at  the  waterline  using  RF  conducting  tape.  This  forms 
a target-image  pair.  This  configuration  was  then  placed  in  a string 
holder  used  for  free  space  swept  frequency  measurements.  A bi-static 
angle  of  approximately  35°  was  used.  Magnitude  and  phase  information 
were  recorded  covering  the  2150  to  4000  MHz  frequency  range  for  both 
vertical  and  horizontal  polarization  (results  shown  in  Figures  58  and 
59  respectively).  Polarization  is  with  respect  to  the  plane  of  the 
target-image  junction.  The  bistatic  separation  is  in  the  plane 
orthogonal  to  the  polarization  vectors.  For  horizontal  polarization, 
three  regions  (large  response  resonances)  are  of  interest  for  both 
bow-on  and  stern-on  incidence  while  nothing  really  significant  is  seen 
at  broadside.  Virtually  the  same  situation  exists  for  vertical  polari- 
zation. It  is  true  that  certain  modes  are  missing  with  this  type  of 
modeling  but  these  would  be  included  with  the  target  on  the  ground  plane. 
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TABLE  II 

COMPLEX  POLES  FROM  SHIP  DATA* 

(Listed  in  order  of  decreasing  residue  magnitude) 


Ship 

Bow 

Stern 

Starboard 

Midway  (1/500) 

-5.83ij26.1 

8.01 ±j75.6 
-147ij0.00 

-18.2ij34.9 
-36.7tj81 .5 
-5.52±jl26 
26.2±j29.2 

-2.13ij27.2 
-1 . lOijO.OO 
-2.42ij74.0 

Missouri  (1/500) 

-7.93ij9.70 

-11.7ijl07 

-5.8Uj5.98 

2.89±jl35 

0.795±jl3.2 

-82.7ij0.00 

-15.4±j87.4 

-3.65ijl2.4 

-124ij0.00 

-11.3ij79.1 

Sverdlov  (1/500) 

-7.3Uj26.0 

5.84ij0.00 

-1.13±jl01 

-0.168±jl3.5 

-12.9±j71.7 

-66.4±jl92 

-26.4ijl8.2 

2.09ij0.00 

-12.9ij95.1 

Missouri  (1/700) 

-0.81Uj8.78 
-1.83±j42.3 
-28.01] 137 

-1 .05tj9.29 
-2.02ij0.00 
-4.99ij96.5 

1 

-19.1ij0.00 

-50.8ij0.00 

-18.3ij74.1 

0.320ij0.00 

Bismark  (1/700) 

-15.6±j27.0 

-24.2ij98.9 

3.15±j74.7 

15.5±jl8.0 

-16.4ij0.00 
-0.753ij33.3 
-30.2ij89.8 
: 7.24ij0.00 

-3.43ijl8.6 

-63.4tj0.00 

-24.2ij61.7 

Mogami  (1/700) 

-5.18ijl5.6 

-146-jl37 

-4.46±j28.1 

-0.938ij7.03 

-37.0ij0.00 

-2.89ij88.6 

-4.90±j8.33 

-3.38±j45.2 

-31.3-jl37 

Hayanami  (1/700) 

-1.63±jl9.3 

-2.15±j38.9 

-3.27±j97.7 

-1.66ijl9.1 

-1.20ij41.5 

-5.37tj91.0 

-10.3ij39.5 

-1.88±jl7.5 

-10.5±j65.2 

-2.84ij97.6 

Shimokaze  (1/700) 

-5.76ij28.8 

-4.97ij50.9 

-0.724±j95.6 

2.92ij0.00 

5.85ij89.3 

-6.33ij30.6 

-44.4ij40.8 

-8.42ij95.1 

-16.7ij20.6 

-25.8ij65.4 

-1.60ij54.8 

-2.96ij99.2 

*The  poles  are  scaled  such  that  the 

imaginary  parts  are 

in  MHz. 
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Figure  50(a) 


76 


DflTfl  FILE 
DflTfl  ID 
SHIFT 
NSTflRT 
NLflST 
• POLES 
VSCRLE 
NRMS  ERR 


SHP28 

MIDWAY  (1/50D) -STERN 
0.0 
1 

200 

8 

11.751 

0.029719 


/\  /M 


200  300  f 

SAMPLE  \ 


Figure  50(b) 
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Figure  52(b) 
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Figure  53(a) 
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Figure  53(b) 
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Figure  53(c) 
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Figure  54(b) 
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Figure  55(b) 


D 


DOTn  FILE!  SHP6 

ORTfl  ID  t HOGRHI (1/700) -STRRBQRRO 
L SHIFT  t 0.0 

NSTflRT  I 1 

I NLRST  I 200 

^ • POLES  I 5 

VSCRLE  I 33.1171 

i MRM<;  FRR  I n.QSQRRR 


OflTfl  FILEi  SHP9 

ORTR  ID  t HRYRNRHI (1/700) -STARBOARD 

SHIFT  I 0.0 

NSTRRT  I 1 

NLRST  t 200 

• POLES  I 8 

VSCRLE  t 14.073 

NRNS  ERR  t 0.002537 


[ 


i 


1 


DATA  FILE 
DRTR  ID 
SHIFT 
NSTRRT 
NLRST 
• POLES 
VSCRLE 
NRMS  ERR 


SHP13 

SHINOKRZE  (1/700) -STRRBORRO 
0.0 
1 

200 

8 

7.932 

0.004747 


I 


97 


OflTfl  FILE 
DATA  10 
SHIFT 
NSTART 
NLAST 
• POLES 
VSCHLE 
NRHS  ERR 


SHP12 

SHINOKAZE  (1/700)  -STERN 
0.0 
1 

200 

6 

0.376 

0.017830 


Figure  57(b) 
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Figure  57(c) 
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gure  59.  Missouri  with  image  (2150-4000  MHz).  Vertical  polarization. 


At  this  time  it  is  planned  to  make  such  measurements  on  the  ground  plane 
range  used  for  the  multifrequency  measurements.  Only  slight  modifi- 
cations are  anticipated.  These  measurements  should  point  to  promising 
frequency  regions  where  the  resonance  behavior  can  be  related  to  sub- 
structures on  the  targets. 

VII.  IDENTIFICATION  OF  NAVAL  VESSELS 

In  this  section  results  of  predictor-correlator  tests  for  dis- 
crimination of  the  naval  vessels  given  in  Table  I are  presented.  The 
signaling  waveform  tested  is  the  matched  filter-type  response  presented 
in  Section  VI.  Note  that  with  this  approach  the  actual  discrimination 
is  based  on  amplitude  data  alone  and  because  actual  measured  data  are 
used  noise  is  present  in  the  signals.  For  each  target  the  classification 
parameters  are  the  complex  poles  extracted  from  the  matched  filter 
responses.  Thus  with  the  present  data  set  each  target  has  four  sets  of 
parameters  corresponding  to  the  aspects  tested  (except  the  Bismark  and 
Mogami  for  which  port  view  data  were  not  taken).  Admittedly  this 
approach  does  not  precisely  satisfy  an  excitation  invariant  criterion. 

It  is  postulated  however  that  each  set  of  parametes  is  quasi-invariant 
in  that  aspect  angles  in  the  vicinity  of  those  tested  will  also  be 
identifiable  from  the  same  four  sets  of  parameters.  Bounds  on  the 
acceptable  aspect  ranges  must  await  additional  measurements.  The  results 
to  be  given  then  are  in  one  sense  a i.'onfirnation  of  the  waveform 
characterization  noted  earlier.  This  itself  appears  to  be  an  advance 
in  the  state-of-the-art.  Our  objective  however  remains  the  exploitation 
of  the  substructure  resonances  of  each  vessel  and  future  studies  are 
directed  to  that  goal.  The  initial  discrimination  results  presented 
here  are  by  no  means  the  optimum  honing  of  the  complex  natural  resonance 
viewpoint. 

In  previous  reports  on  discrimination  testing  volumous  plots  of 
the  correlation  as  a function  of  the  sample  density  (Equation  (13))  have 
been  given,  A more  compact  presentation  of  such  results  is  shown  here. 
Each  harmonic  data  set  for  a given  vessel  and  orientation  was  tested 
(Equation  (13))  against  all  vessels  and  aspects  via  the  complex  poles 
in  Table  II  (data  for  all  four  aspects  for  each  vessel  were  used). 

From  each  of  these  computations  two  numbers  were  extracted.  These  were 
the  average  value  of  the  correlation  coefficient  and  the  minimum  value 
of  the  correlation  coefficient.  The  average  value  is  defined  by 

N 

^ p"(nAt) 

<P’’>  = (46) 
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where  the  sumnation  is  over  all  sample  densities.  A search  was  done  for 
each  of  these  correlation  parameters  for  each  vessel  aspect  to  obtain 
the  pole  set  which  yielded  the  largest  value  for  each  parameter.  Thus  two 
simple  identification  parameters  were  tested.  Two  example  plots  are 
shown  in  Figures  60  and  61.  In  both  cases  the  harmonic  data  for  the 
Missouri  (1/500)  as  viewed  from  bow-on  was  used.  In  Figure  60  the 
poles  of  the  above  data  set  were  used.  In  Figure  61  the  poles  taken 
from  the  Shimokaze  as  viev/ed  from  the  stern  were  used.  Table  III 
summarizes  the  results  based  on  predictor-correlator  processing  of  thirty 
harmonic  data  sets  and  thirty  complex  pole  sets  of  all  aspects  of  the 
ships  in  Table  I.* 


TABLE  III 


■ - - 1 

Method 

Probability  of  Correct 
Classification 

Maximize  <p"> 

0.77 

Maximize  (p"^„) 
mm 

0.70 

There  is  obviously  some  risk  involved  in  results  based  on  a single 
number  as  opposed  to  thresholds  selected  over  a range  of  sample 
densities.  Nevertheless  Table  III  is  felt  to  be  a fair  estimate  of  our 
present  capability.  Figure  62  shows  the  results  of  the  predictor- 
correlator  processing  of  all  pole  sets  with  the  harmonic  data  of  the 
(1/500)  scale  Missouri  viewed  from  bow  on.  The  solid  curve  (actually 
the  marked  points)  is  for  <p">  and  the  dashed  line  for  (prnin)-  For 
these  harmonic  data  the  correct  pole  set  was  number  one  and  it  can  be 
seen  that  both  parameters  yielded  correct  identification.  It  might 
appear  that  (pmin)  would  be  a more  effective  discrimination  tool  based 
on  this  isolated  example.  This  may  be  somewhat  misleading  since  based 
on  the  results  of  Table  III  it  is  felt  that  <p">  is  the  more  effective 
parameter  of  the  two.  An  example  of  this  is  shown  in  Figure  63  where 
the  harmonic  data  for  the  Mogami  as  viewed  from  the  stern  has  been 
compared  to  all  pole  sets.  For  these  data  the  correct  pole-set  is 
number  19.  For  this  example,  identification  based  on  (Pnip)  fails 
since  pole-set  number  16  was  chosen.  Identification  based  on  <p">, 
however  was  successful.  It  is  also  felt  that  the  averaging  procedure 
should  be  less  susceptible  to  noise  than  the  one-point  minimization 
procedure. 


*Tlie  thirty  data  sets  represent  four  views  of  all  ships  in  Table  I 
except  the  Bismark  and  the  Mogami  for  which  port  data  was  not  taken. 
The  poles  taken  from  the  port  view  data  are  not  listed  in  Table  II 
since  they  are  nearly  identical  to  those  for  the  starboard  views. 
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Missouri  (1/500)  - Bow  On  Harmonic  Data 

Shimokaze  (1/700)  - Stern  Poles 
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VIII.  CONCLUSIONS  AND  RECOMMENDATIONS  j 

A.  Conclusions 

1.  A method  for  discrimination  of  naval  vessels  using  | 

amplitude  scattering  data  at  ten  discrete  iiamonical ly 

related  frequencies  has  been  described  and  illustrated 
using  measured  radar  data  tests  over  five  vessel  classes 
and  four  widely  separated  aspects  for  each  vessel  yields 
a 77%  probability  of  identification. 

2.  The  discrimination  achieved  in  conclusion  1)  uses  j 

only  amplitude  data  and  noise  is  inherently  present  in 

the  measurements.  The  method  is  far  from  optimum  and 

does  not  completely  exploit  an  excitation  invariant  set  j 

of  target  descriptors.  The  initial  results  therefore,  1 

while  an  advance  in  vessel  identification,  are  merely 

indicative  of  what  might  be  achieved  with  precise,  I 

excitation  invarient  substructure  resonances.  1 

3.  A ten  frequency  reflectivity  range  with  the  sea  medium 
modeled  as  a perfectly  conducting  ground  plane  has  been 
instrumented  and  data  collected.  Some  nine  models  of  various 
vessel  classes  have  been  assembled  and  measured  at 
selected  aspects. 

4.  A general  philosophy  for  the  identification  of  naval 
vessels  based  on  the  complex  natural  resonances  of 
selected  substructures  of  the  vessels  has  been  described. 

Special  attention  has  been  given  to  the  inherent  problems 
of  naval  vessels;  large  electrical  size  and  the  presence 
of’an  interfering  sea  medium. 

5.  Studies  of  methods  for  extracting  desired  complex 

natural  resonances  from  measured  scattering  data  have  [ 

been  initiated  and  some  progress  achieved.  In  the  presence  I 

of  noise  and  with  weakly  excited  resonances  the  problem 

is  difficult  but  not  insurmountable.  A real  advantage  of  r 

the  techniques  described  in  this  report  is  that  the  j 

resonances  need  not  be  extracted  in  real  time.  That  is,  ‘ 

the  discrimination  procedure  does  not  extract  resonances 

from  the  signaling  waveform  of  an  unVnown  target.  j 
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IX.  RECOMMENDATIONS 

1.  Swept,  in  addition  to  discrete,  frequency  measurements 
appear  desirable  for  naval  vessel  targets.  With  swept  frequency 
measurements,  synthetic  "chirp"-type  processing  will  permit 
isolation  of  potentially  interesting  structure  features.  At  the 
same  time  the  measured  data  will  suggest  possible  high  return 
frequency  regions.  Modifications  to  convert  the  existing 
reflectivity  facility  to  discrete  or  swept  measurements  are  in 
progress. 

2.  Improvements  are  clearly  still  needed  in  methods  for  extracting 
the  complex  natural  resonances  of  a target  from  measured  response 
data.  It  is  not  intended  that  the  main  purpose  of  this  program  be 
directed  to  this  goal  (the  same  question  arises  on  other  programs) 
but  some  study  will  be  made  of  alternative  schemes,  one  such 
being  via  linear  algebra  to  find  that  difference  equation  which 
satisfies  several  data  sets. 

3.  A method  for  obtaining  many  complex  natural  resonances  from 
measured  data  obtained  with  an  internally  located  (within  the 
target)  source  and  receiver  has  been  suggested  [18,19]. 

Experimental  tests  of  this  procedure  will  be  conducted  during 
the  next  interim. 

4.  Additional  measured  scattering  data  on  model  naval  vessels, 
in  situ,  are  needed.  These  include  swept  frequency  measurements 
at  untested  aspects  and  possibly  additional  classes  of  naval 
vessels. 

5.  Assuming  that  research  progress  on  the  identification  of  naval 
vessels  continues  to  show  workable  methods,  some  study  should  be 
initiated  on  the  actual  size  and  complexity  of  full  scale  radar 
systems  needed  for  the  identification  task.  Ideally  such  a study 
should  be  started  late  in  the  upcoming  interim. 
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